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ABSTRACT 

Context. The INTEGRAL observatory operating in a hard X-ray/gamma domain has gathered a large observational data set over nine 
years starting in 2003. Most of the observing time was dedicated to the Galactic source population study, making possible the deepest 
Galactic survey in hard X-rays ever compiled. 

Aims. We aim to perform a Galactic survey that can be used as the basis of Galactic source population studies, and perform mapping 
of the Milky Way in hard X-rays over the maximum exposure available at \b\ < 17.5°. 

Methods. We used sky reconstruction algorithms especially developed for the high quality imaging of INTEGRAL/IBIS data. 
Results. We present sky images, sensitivity maps, and catalogs of detected sources in the three energy bands 17 — 60, 17 — 35, and 
35 — 80 keV in the Galactic plane at \b\ < 17.5°. The total number of sources in the reference 17 — 60 keV band includes 402 objects 
exceeding a 4.7cr detection threshold on the nine-year time-averaged map. Among the identified sources with known and tentatively 
identified natures, 253 are Galactic objects (108 low-mass X-ray binaries, 82 high-mass X-ray binaries, 36 cataclysmic variables, 
and 27 are of other types), and 115 are extragalactic objects, including 112 active galactic nuclei (AGNs) and 3 galaxy clusters. The 
sample of Galactic sources with S /N > A.lcr has an identification completeness of ~ 92%, which is valuable for population studies. 
Since the survey is based on the nine-year sky maps, it is optimized for persistent sources and may be biased against finding transients. 
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1. Introduction 

A large fraction of astrophysical phenomena cannot be studied 
via observations of individual sources, but require instead large 
statistical studies. The last few decades have provided us with 
great opportunities for studies of the populations of compact 
sources (black holes, neutron stars, white dwarfs) in our Galaxy 
and nearby galaxies. 

In particular, surveys of the sky in hard X-rays were per- 
formed with the IBIS telescope ( Ubertini et al. 2003 1 of the 
INTEGRAL obs ervatory (, Winkler et al.| |2003| l and Burst Alert 



Telescope (BAT; Barthekny et al. 200511 at the Swift observa 



tory (Gehrels et al. 2004 1. In contrast to Swift, with a nearly 



uniform all-sky survey, which is especially useful for studies of 
active galactic nuclei (AGN; [Tueller et al.| |2010[ |Cusumano et 
[aL] [20Tbl [Ajello et al.| [20T2| , the INTEGRAL observatory pro- 
vides a sky survey with exposure that are deeper in the Galactic 
plane (GP) and has higher angular resolution, which is essen- 
tial in these crowded regions. This makes the Swiftj'QPsJ and 
INTEGRAL/IBIS surveys complementary to each other. 

The INTEGRAL observatory has been successfully operat- 
ing in orbit since its launch in 2002. Over the past few years, 
INTEGRAL data has allowed us to construct high quality cat- 
alogs (Revnivtsev^eTaLl |2004b| [2006) [Molkov et al.[ [2004 



IKrivonos et al.l |20051 |2007b| |2010bt [Bird et al.l |2006| 12007 



* Based on observations with INTEGRAL, an ESA project with in- 
struments and science data centre funded by ESA member states (espe- 
cially the PI countries: Denmark, France, Germany, Italy, Switzerland, 
Spain), Czech Republic, and Poland, and with the participation of 
Russia and the USA 



2010[l, to reveal n ew types of sources ( Courvoisier et al. 2003 
Revnivtsev et al. 2003d Sguera et al. 2006)1, t o calculate the 



statistics of active galactic nuclei ( Sazonov et al. 

|2009f : 



20071 Bassani 



|et al.| |200 6; Beckma nn et al.| |2009| ), of low mass X-ray bi- 
naries ^Revnivtsev et al. , 2008a), of high mass X-ray binaries 



( [Lutovinov et al.| |2005b, ,2007, |Bodaghee et al.j [20071 |20T2ll, 
and of catac lysmic variables ( [Revnivtsev et al.| |2()08b[ |Scaringi| 
^taLl[20T0l i. 



In our previous papers (Krivonos et al. 2010a|b I, we pre- 
sented the seven-year hard X-ray all-sky survey in the energy 
range 17 — 60 keV based on the improved sky reconstruction 
method for the IBIS telescope. The sensitivity of the survey was 
significantly improved by suppressing the systematic noise. 

Here we present the selected nine-year averaged sky im- 
ages, sensitivity maps, and catalog of the sources detected in the 
Galactic plane {\b\ < 17.5°) in three energy bands: 17 — 60, 
17 — 35, and 35 — 80 keV. This survey is the most sensitive X- 
ray survey of the Galaxy existing so far and it can be used as: 1) 
a basis for statistical studies of different types of sources in our 
Galaxy, and 2) a guiding line for new surveys of a new genera- 
tion of hard X-ray focusing telescopes (e.g. NuSTAR described 
in [Harrison etaL[ ( [20T0[ ), and Astro-H in |Takahashi et aI.| ( [20T0| l). 

The full se t of sky maps is available at the Sky View Virtual 
Observator}[^(McGIynn et al. 1998 1 and Russian Science Data 
Cente:[^for the INTEGRAL observatory at the Space Research 
Institute (IKI), Moscow. 



http://skyview.gsfc.nasa.gov 
http://hea.iki.rssi.ru/integral 



1 



Krivonos et al.,: INTEGRAL 9-year Galactic Survey 




Fig. 1. Dead time-corrected exposure map of the INTEGRAL 
all-sky survey (January 2011, public data). The grid gap in the 
Galactic latitude is 17.5°, which is a half-height of the current 
Galactic survey. Blue contours represent exposure levels of 10, 
150, 800, 2000, and 4000 ks. The effective exposure in the GC 
region is 12 Ms, which corresponds to 26 Ms of a nominal time. 



2. Survey 

To conduct the current Galactic survey, we selected pub- 
licly available INTEGRAL data from December 2002 to 
January 2011 (spacecraft revolutions 26-1013). Every individ- 
ual INTEGRAL observation with typical exposure time of 2 ks 
(so called Science Window, ScW) was analyzed with a specially 
developed software package (see e.g. Kri vonos et al.|2010a| and 
references therein) to produce sky images in three energy bands: 
17 — 60, 17 — 35, and 35 — 80 keV. In contrast to our previ- 
ous surveys, the flux scale in each ScW sky image was adjusted 
using the flux of the Crab nebula measured in the nearest ob- 
servations. This procedure was used to account for the ongoing 
detector degradation and loss of sensitivity at low energies. 

In total, we obtained 73489 sky images in each band, which 
comprises ~ 132 Ms of the effective (dead time-corrected) ex- 
posure. The survey sky mapping was organized in six over- 
lapping 70° X 35° Galactic cartesian projections centered on 
zero Galactic latitude (b = 0°) and / = 0°, +50°, + 115°, and 
I = 180°. The latitude coverage of the survey \b\ < 17.5° was 
chosen with the IBIS/ISGRI field of view (28° x 28°) and stan- 
dard 5° X 5° observational pattern in mind. Thus, we used all 
observations performed by INTEGRAL in the Galactic plane. 
Fig.[T]illustrates the INTEGRAL exposure map of publicly avail- 
able observations up to January 201 1. 

The survey sky coverage versus the 4.7cr limiting flux is 
shown in Fig. [2] The peak sensitivity of the survey is 2.9 x 
10"'^ erg s"' cm~2 (- 0.20 mCrab in 17-60 keV) at a 
4.7cr detection level. The survey covers 90% of the geomet- 
rical area (12680 degrees) down to the flux limit of 2.0 x 
10"" erg s"' cm-2 (~ 1.41 mCrab) and 10% of the total 



10 erg s ' cm ^ 



area down to the flux limit of 4.9 x 
0.30 mCrab). 

2.1. Systematic noise 

INTEGRAL/IBIS deep sky mosaics are usually affected by a 
systematic noise caused by the source confusion in the region 
of GC and by the imperfect sky reconstruction ( [Krivonos et al.| 
|2010a| l. 

The large field of view (FOV) of the INTEGRAL/IBIS tele- 
scope leads to a high probability of having some bright X- 
ray source within the instrument FOV during any galactic ob- 
servations. Therefore, to work at the level of Poisson noise, 
the INTEGRAL/IBIS telescope image reconstruction procedure 




10-" lO-'o 
4.7(7 Limiting flux. 17 — 60 keV (erg s"' cm"^) 

Fig. 2. Sky area as a function of the limiting flux for source de- 
tection with a 4.7cr significance (solid line). The black dashed 
curve shows the sky coverage with the masked area around 
bright sources listed in Table [T] The dotted line represents the 
geometrical area. 

Table 1. Exclusion radius around bright sources. 



Name 


Exclusion radius, deg. 


Crab 


23.2 


Sco X-1 


14.0 


Cyg X-1 


4.8 


Cyg X-3 


3.9 


Vela X-1 


3.3 


GX 301-2 


2.2 


GRS 1915+105 


3.2 



should have a dynamic range of 10^ and more, which is very dif- 
ficult to achieve owing to the imperfect modelling of the mask 
shadow, the individual pixel sensitivity, the variability of the 
background pattern, etc. 

In spite of the latest version of the IBIS sky image recon- 
struction allowing us to reach the dynamic range of the images 
~ 10^, some sky artefacts are still present around bright sources, 
such as Crab, Sco X-1, Cyg X-1, Cyg X-3, Vela X-1, GX 301-2, 
and GRS 1915+105 (see a sky mosaic at / = +50° in Fig.|3j. 
To prevent any false detections, we masked out circular regions 
around bright sources with radii listed in Table [T] The exclu- 
sion radius was chosen to contain all significant (> 4.7<t) neg- 
ative excesses (indicators of the systematic noise, which is as- 
sumed to be symmetric) around a given bright source. We note 
that known sources with S /N > lOo" falling inside these re- 
gions were included in the source catalog. A rejection of areas 
with high systematic noise significantly improved the quality of 
the survey mosaics, which is demonstrated in Fig. |4] where we 
show a distribution of signal-to-noise (S /N) values for pixels in 
the sky mosaic at / = +50°. As seen from Fig.|4] the masked 
sky image does not contain strong deviations from the Gaussian 
distribution, in contrast to the original one. Masked areas around 
bright sources reduce the geometrical area of the survey by 13% 
(Fig.|2). 

Systematic effects are less important in the harder energy 
bands. Fig. |5] shows that a sky image of the GC region in the 
35 — 80 keV energy range is practically free of the systematic 
noise with respect to 17 — 60 keV band, which is confirmed by 
the S /N distribution of image pixels in Fig.|6] 
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55.0' 50.-D' 45.0' 

Galactic Litngltude 



Fig. 3. INTEGRAL/IBIS hard X-ray (17-60 keV) map of the sky region of Cyg X-1, Cyg X-3, and GRS 1915-)-105 at I 
with masked area shown as green circles. The corresponding S /N distribution of pixels is shown in Fig.|4] 



+50° 



100 




gion contains 12 peaks above 4.7cr in the reference 17 — 60 keV 
energy band. All these candidate sources have a very distorted 
shape, and none of them have been detected in the 35 — 80 keV 
energy range, which also points to a false detection. Therefore, 
all excesses in the GC region have been attributed to the system- 
atic noise. 



100 



■10 10 

Signal-lo-noise ratio 



Fig. 4. Signal-to-noise ratio distribution of a number of pixels in 
the 17 — 60 keV map of Fig. [3] Red and blue histograms show 
distributions of pixels in full and masked images, respectively. 
The dashed line represents the normal distribution with unit vari- 
ance and zero mean. 



2.2. Detection threshold 

Given the IBIS/ISGRI angular resolution of the survey, the sky 
map contains ~ 3 x 10^ statistically independent pixels. Taking 
into account the minor contribution of the systematic noise, we 
adopted a conservative detection threshold of {S /N)\im > ^-.Tcr 
to ensure that the final catalog contains no more than one spuri- 
ous source assuming Poisson statistics. 

A source detection in the region of ~17 degrees around 
GC should be interpreted with special care because of the pos- 
sible false peaks induced by the systematic noise. The latter 
is revealed by a distorted shape that differs significantly from 
the instrumental point-spread function (PSF), which is a two- 
dimensional Gaussian (cr = 5')- The sky map of the GC re- 




-10 -5 5 

Signal-lo-noise ratio 

Fig. 6. Signal-to-noise ratio distribution of a number of pixels in 
the hard X-ray image shown in Fig.|5] The dashed line represents 
the normal distribution with unit variance and zero mean. Red 
and blue histograms show pixel distributions for images in the 
energy band 17 — 60 keV (high systematic noise, see text) and 
35 — 80 keV, respectively. 



3. Catalog of sources 

The catalog was compiled from the source sample exceeding the 
detection threshold in the reference 17 — 60 keV energy band. 
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Fig. 5. INTEGRAL/IBIS hard X-ray (35 - 80 keV) map of the sky region around the GC. The total exposure is about 26 Ms in the 
GC region. The image is shown in terms of S /N with the color map ranging from values of to 25. This color scheme is used to 
emphasize sky background variations. Fig. [6] demonstrates a corresponding S /N distribution of pixels. 



The list of sources is presented in Table 2, and its content is 
described below. 

Column (1) "Id" - source sequence number in the catalog. 

Column (2) "Name" - source name. Their common names 
are given for sources whose nature was known before their de- 
tection by INTEGRAL. Sources discovered by INTEGRAL or 
those whose nature was established thanks to INTEGRAL ob- 
servations are named "IGR" 

Columns (3,4) "RA, Dec" - source equatorial (J2000) coor- 
dinates. 

Column (5,6,7) "Flux" - time-averaged source flux in the 
17-60 keV, 17-35 keV, and 35 - 80 keV energy bands, re- 
spectively. 

Column (8) "Type" - general astrophysical type of the ob- 
ject: LMXB (HMXB) - low- (high-) mass X-ray binary, AGN 
- active galactic nucleus, SNR - supernova remnant, CV - cat- 
aclysmic variable, PSR - isolated pulsar or pulsar wind nebula 
(PWN), SGR - soft gamma repeater, RS CVn - coronally ac- 
tive binary star, SymbStar - symbiotic star. Cluster - cluster of 
galaxies. A question mark indicates that the specified type is not 
firmly determined, so should be confirmed. 

Determining of the natures of the sources is complicated and 
a continually ongoing process. During the compilation of the cat- 
alog, we selected from our point of view, the most recent or re- 
liable, identifications from the literature. In some cases (when 
references are not given), the identifications had been performed 
using recent observations of the Chandra, XMM-Newton, and 
Swift observatories, both the Simbacj^and NEE(^ database, as 
well as the 2MASS catalog. 

Column (9) "Ref." - references. These are mainly provided 
for new sources and are related to their discovery and/or nature. 
The papers describing a routine analysis of a given source e.g. 

^ Simbad Astronomical Database http://simbad.u-strasbg.fr 

* NASA/IPAC Extragalactic Database http://ned.ipac.caltech.edu 



confirmation of their nature, a refined position, etc. are not ref- 
erenced. 

Column (10) "Notes" - additional notes such as type sub- 
class, redshift information, alternative source names. The red- 
shifts of the extragalactic sources were obtained from the 
Simbad and NED databases. 



4 



Krivonos et al.,: INTEGRAL 9-year Galactic Survey 



4. Concluding remarks 

Our Galactic survey is based on sky maps averaged over a nine- 
year time period, which obviously means that it has a strong bias 
against finding low S /N transient sources. This ensures that the 
current survey contains mainly persistent objects, in addition to, 
however, objects with strong intrinsic variability. 

We have presented detailed source statistics in Table [3] and a 
chart of source types in Fig.|7] Our Galactic survey has an iden- 
tification completeness of {Njot—NNotiDj/NTot = 1—34/402 
0.92, which provides a strong statistical basis for population 
studies. In the complementary paper of Lutovinov et al., (2012, 
in prep.), we note in particular that we used the source sample 
and sensitivity maps of the current survey to study the HMXB 
luminosity function and their spatial distribution in the Galactic 
disk. 



o 




AGN LMXB HMXB CV Other NotID 



Fig. 7. Chart for the source classes detected at S /N > 4.1 cr in 
the reference 17 — 60 keV energy range (see Table[3]). Green bar 
denotes the number of sources that have a tentative association 
with a given type of objects. 
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Table 3. Census of Galactic and extragalactic sources at \b\ < 17.5° 



Energy range, keV 


AGN 


LMXB 


HMXB 


CV 


Other 


NotID 


Total 


17-60 
17-35 
35 - 80 


104+8* 

89+7* 

76+4* 


100+8* 

98+7* 

80+5* 


76+6* 
74+6* 
66+3* 


35+1* 
35+1* 
13 


30 
28 
26 


34 
22 
7 


402 
367 
280 





Notes. The number of sources with a tentative classification is denoted with S index. Fig. |7] shows the chart for source classes detected in the 
17 — 60 keV energy band. 
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